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■ Abstract. We calculate some 0(a 2 s ) corrections to the JIMWLK kernel in the 
framework of the light-cone wave function approach to the high energy limit of QCD. 

f- ^ The contributions that we consider originate from higher order corrections in the strong 

[ — . coupling and in the density of the projectile to the solution of the classical Yang-Mills 

equations of motion that determine the Weizsacker- Williams fields of the projectile. 

■ We study the structure of these corrections in the dipole limit, showing that they are 
subleading in the limit of large number of colours N , and that they cannot be fully 

O i' recast in the form of dipole degrees of freedom. 

<d : 

^ . 1. Introduction: the wave function formalism and JIMWLK 

S : 

The study of high gluon density QCD, pioneered in the MV model PQ, has become a 
fashionable subject both due to its theoretical relevance and to eventual applications to 
small-x DIS or high-energy nuclear collisions (see the review [2]). The main outcome 
of these studies was the JIMWLK equation which controls the evolution of the colour 
distribution of a hadron with increasing energy. In the last 3 years, the limitations 
of such evolution equation have become apparent in the form of the discussion of 
fluctuations and pomeron loops, see the reviews [3J. In all available schemes, a separation 
of scales is considered between fast partons which act as sources for the classical 
dynamics of the soft, slow glue. Considering gluon radiation and the arbitrariness of the 
fast-slow separation led to a renormalization group equation, the JIMWLK equation. A 
powerful framework to derive this equation, in which its limitations become apparent, 
is the wave function approach (WFA) [4J which we briefly introduce in this Section, see 
full details in [HE]. 

In the WFA the evolution is introduced in the wave function of the projectile 
scattering off a hadronic target. This wave function is the superposition of fast gluons 
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with k + > A. The scattering matrix with the target is the superposition of independent 
scattering matrices S ab . Then, an average over target configurations is taken, 

E Y [S] = (¥(Y)\&(Y)), (£) = J dSZ[S]W Y [S], (1) 

with Wy [S] the probability density for the target to have a certain configuration of the 
fields. After a boost Y — > Y + 8Y, the wave functions become 

\<t'(Y + 6Y)) = {l - \SY J d 2 zbUz,Sp)bUz,Sp) (2) 

with S ab — > 5 ab for |\I/ l (Y + 5Y)), p the projectile density operator and 6" the 
Weizsacker- Williams fields of the projectile (WW) determined from the classical Yang- 
Mills equations of motion (EOMs), in which the hard gluons enter as an external source. 
The resulting evolution equation is JEfSj/^Y = xE[5*], with the kernel given by 
1 
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WlMJl) + b?(J R mjR) - 2S ba {z)b1{J L )b\{J 1 



(3) 



where b\ depends on z and Jl{b) are left (right) colour rotation operators, whose action 
in the dipole model limit, with the projectile composed of colour dipoles [1], is illustrated 
in Figured! With the EOMs solved at lowest order in gp in the A + = gauge p], 

bUz) = 9 (b^z) + gHUz) + 0(g 4 )) , b^z) = j^L ftX) f(z),(4) 
with diX = df In (\z — x\X), the JIMWLK evolution equation is obtained [1]. 
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Figure 1. Diagrammatic representation of J£E (top) and J^J^T, (bottom). The 
target is represented by the vertical thick line. 



2. First correction from the classical equations of motion 

Two limitations of the derivation leading to JIMWLK are apparent: First, Equation (j2j) 
is an approximation only valid for a low-density projectile. Second, the solution of the 
EOMs is done at the lowest order in the density of the projectile, Equation PJ. In this 
Section we explore the corrections to the latter by going to order g 3 p 2 [HIE], 

hUz) = ~4(hT 2fabC J d2xd2 y( X 9iY) P b (x)p c (y), (5) 
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X diY = X (d?Y) — (dfX) Y. With this expression the expansion of the kernel becomes 
X = a.(x 1) + a a x 2) + 0(a%j), (6) 
whose leading term is the JIMWLK kernel 
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mmv) + r R .r R - 2S ba r L j R 



2,7T Jxyz 

with the J's and S depending on z, and the first correction reads 

/2) = -7A2 / f abc (d l Y)(Xd l W) [ji(x)Jt(y)Ji(w) + BBR 
2f acd (d i X)(Yd i W)S ba (z) \jl(x)J R (y)jf t (w) + LLR 



(7) 



X 



(8) 












d 
a 

















1 ^ : 


1 








1 ^ : ; 1 




1 








1 


1 











a 


b 










U '-^ a 




e 




















1 ^ : 


i 



Figure 2. Diagrams corresponding to LLL (plots on the left) and LLR (plots on the 
right), see the text. The diagrams include one (oc 1/N, top), two (oc 1/N 2 , middle) 
and three (oc 1/N 3 , bottom) active dipoles. 



To proceed further we go to the dipole model limit, in which the action of the 
different terms can be explicitly computed and is illustrated in Figure [2j With s(x, y) 
the .S-Hiatrix for a dipole, the kernel (jHJ) can be separated [5] into three pieces containing 
one, two or three active dipoles with the corresponding 1/N suppression factors, 

x 2) ns] = ( X f /N + X ? /N * + X? /N s) £[*]■ (9) 

Xi/n * s found to vanish. Thus the solution to next-to-leading order in the EOMs does 



not yield leading 1/N corrections and, therefore, it does not provide a correction to the 
BK equation [7]. Xi/n 2 anc ^ Xi/n 3 cann °t be fully recast in terms of dipoles. Noting 
that any wave function or weight functional of a gluonic/dipole configuration has to be 
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completely symmetric under the exchange of any number of gluons/dipoles, we get 

X f /N ^[s] = i J dh^R-P^Ud.V) (10) 

-d(V-U){PdiR)][Nix{SlS v S$S r SlS g }]+di(V-U) 

x (R-P) g (V-U) [tr{SiS v SiS r }tr{StS z }-tr{StS z SlS r }tr{StS v } 

~\~to{SpS z Sl l/ S r }t'i:{Sj,S v } — tr {SpS v SlS r }tr{Sl l S z } j — -. 



5s (p, r)5s(ix, u) 

X^aSH = -\j^yJpl d2z [{^(R-P)}{{T~Q) 5 (V- £/")}" 
x |^tr{5'^5' 1 ,5' ( |S't5'pS' r } + tr{S'^S'(S'jS' r .*S'pS , t; } + trlS'^SV^'pS^S'^S^S'gSt} 
+tr{5{5 r 5l5 f ,fllS t 5j5 J( }] ^^^^) ■ ( n ) 

3. Conclusions 

In this contribution we present the 0(a 2 ) corrections to the JIMWLK evolution equation 
coming from the 0(g 3 p 2 ) solution to the classical EOMs. The leading 1/N piece in the 
dipole model limit is found to vanish, thus yielding no correction to the BK equation. 
Subleading corrections do not show a closed dipole form. While the corrections we 
compute are certainly not the complete set of 0{a 2 s ) corrections to JIMWLK, they are 
part of the full solution and may turn to be important to fulfill general requirements of 
the complete theory of high gluon density QCD like dense-dilute duality [H]. 
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